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Abstract 

We present a Q CD-based approach to the endpoint region of the 
lepton spectrum in 5 — > Xu i v decays. We introduce a genuinely 
nonperturbative form factor, the shape function, which describes the 
fall-off of the spectrum close to the endpoint. The moments of this 
function are related to forward scattering matrix elements of local, 
higher-dimension operators. We find that nonperturbative effects are 
dominant over a finite region in the lepton energy spectrum, the width 
of which is related to the kinetic energy of the 6-quark inside the B 
meson. Applications of our method to the extraction of fundamental 
standard model parameters, among them Vub-, are discussed in detail. 



(submitted to Physical Review D) 



1 Introduction 



Recently, much progress has been achieved in the understanding of inclusive 
weak decays of hadrons containing a heavy quark Q. Using the theoretical 
tools of the operator product expansion and the heavy quark effective theory 
(HQET) [|I], H, 1, 1^, 1^ , one can construct a systematic expansion of the (dif- 
ferential) decay distributions in powers of A/mg, where A is a characteristic 
low-energy scale of the strong interactions H, ^ ^ . Quite remarkably, 
the parton model emerges as the leading term in this QCD-based expansion, 
and the nonperturbative corrections to it are suppressed by a factor A^/mg. 
The fact that there are no first-order power corrections relies on a particular 
definition of mg, which is provided in a natural way by requiring that there 
be no residual mass term for the heavy in HQET [|ll|, 12 1. This definition 



is unique and can be regarded as a nonperturbative generalization of the 
concept of a pole mass. 

The availability of a systematic, QCD-based expansion of inclusive decay 
rates raises the hope for a better understanding of these processes in gen- 
eral, and in particular for a more reliable extraction of the standard model 
parameters m^, mc, Kt, and Vuh, which was so far hindered by strong model 
dependence. For a determination of Vub^ however, it is essential to understand 
the endpoint region of the lepton spectrum, which is of genuinely nonper- 
turbative nature. Although the new methods developed in Refs. 0, ||, [T0| 
provide an important step towards this goal, they are not directly applicable 
to the endpoint region. The difficulties arise from the fact that close to the 
endpoint the expansion parameter is no longer A/m;,, but A/(mf, — 2Ei), 
and thus the theoretical prediction becomes singular when the lepton energy 
approaches the parton model endpoint -E^^max = ''^6/2. It is then not obvious 
how to interpret the theoretical predictions. 

To see what the problem is, consider the theoretical prediction for the lep- 
ton spectrum in i? — ^ t v decays. Including the leading nonperturbative 
corrections, one obtains |, 



where 

nil) ' 1927r 



ml , (2) 



1 



and 



F{y) = (3 - 2y) y + ^ + (6 + 5y) y ^ . (3) 
For simplicity, we do not include perturbative QCD corrections, which have 



been calculated in Refs. [|T^, 0. The hadronic parameters Ai and A2 are 
related to the kinetic energy Kh of the heavy quark inside the B meson, and 
to the mass splitting between B and B* mesons [|15|: 

K, = -^, ml,-ml = 4X,. (4) 

The singular structure of the operator product expansion close to the end- 
point at y = 1 manifests itself in the appearance of (5-function (and higher) 
distributions. Certainly, one cannot trust the shape of the theoretical spec- 
trum close to the endpoint. Nevertheless, integrating (P with a smooth 
weight function, one obtains well-behaved results for quantities such as the 
total decay rate and the average lepton energy: 

2ml 20 I Uml ' ^ ^ 



The coefficients of the singular terms give nonvanishing contributions to these 
integrated quantities. There is thus some relevant physical information con- 
tained in these terms. 

Bigi et al. have advocated to integrate over the singularities before 
confronting the theoretical prediction for the lepton spectrum with data. 
They argued that one has to integrate over a finite energy interval of at least 
several hundred MeV, corresponding to a region of order 1 / rrib in the variable 
y. This proposal is based on the idea of quark-hadron duality, which implies 
that when one sums over a sufficient number of exclusive hadronic modes, 
the decay probability into hadrons equals the decay probability into free 
quarks.^ Note that the 5-function term in (|I|) contributes to the integrated 
spectrum, but the term proportional to 6'{1 — y) does not. Similarly, more 



^One expects that duality holds for the lepton spectrum even in the endpoint region, 
which extends over an interval of order 1/mi, in y. Only in a tiny region of order l/?n^ 
below the physical endpoint, the spectrum is dominated by a few exclusive modes. 
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singular terms, which appear at higher orders in the l/nii, expansion, do not 
contribute. 

A shghtly modified procedure was proposed by Manohar and Wise [^j. 
They chose to smear the spectrum with a gaussian distribution of width Ay. 
Empirically, they found that Ay ~ 0.2 — 0.5 is necessary to obtain from the 
theoretical prediction a smooth lepton spectrum, which can be compared to 
data. This procedure has the disadvantage that the results depend on the 
smearing function, and that the choice of Ay is ad hoc. When the smearing 
function is chosen to be symmetric, it again follows that the term proportional 
to 5'{\ — y) does not contribute to the smeared spectrum. 

The frustrating conclusion of these analyses is that the new theoretical 
methods are only of very limited use for a more reliable determination of Vub, 
since the region of the lepton spectrum which is accessible to a measurement 
is smaller than the region over which the theoretical spectrum has to be 
integrated in order to obtain a reasonable result. 

As proposed in Refs. 0, H, the theoretical description is to a large 
extent ignorant to the rich physical information contained in the lepton spec- 
trum close to the endpoint. In this paper, we shall suggest a different ap- 
proach. It is motivated by a very simple observation: In the parton model, 
the endpoint region of the lepton spectrum is described by a step function, 
the location of which is determined by the kinematics of a free quark decay. 
The true physical endpoint, however, is determined by the decay kinemat- 
ics of hadrons. Hence, when QCD is trying to tell us something about the 
redistribution of the endpoint region due to nonperturbative effects, it can 
only do this by the occurance of singular functions. Our approach will allow 
us to extract the physical information contained in the singular terms in the 
QCD-predicted lepton spectrum in a systematic way, and to all orders in 
the l/rrih expansion. To this end, we shall introduce the concept of a shape 
function, which is a genuinely nonperturbative form factor that describes the 
fall-off of the spectrum in the endpoint region. We find that the character- 
istic width of this region is given by ay = (— Ai/3m^)^/^, corresponding to a 
finite region in the lepton energy. Although there do not appear first-order 
power corrections in (|l|), it is important to realize that there exists a small 
region where the true spectrum is very different from the theoretical predic- 
tion. This difference is described by the shape function. We will show that 
the moments of this function can be addressed in QCD, and can be related 
to hadronic parameters (such as Ai and A2) that are defined in terms of for- 
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ward scattering matrix elements of local, higher- dimension operators. To all 
orders in 1/mi,, the leading contributions to the moments can be given in 
closed form. 

We believe that our approach will eventually lead to a better under- 
standing of the nonperturbative aspects of the lepton spectrum close to the 
endpoint. It establishes the connection between the experimentally observed 
lepton spectrum and the underlying theory of QCD. This connection works 
in both directions: Theoretical ideas about the moments of the shape func- 
tion can help to analyze the lepton spectrum and to determine the values of 
the quark masses and mixing angles. On the other hand, from a precise mea- 
surement of the spectrum in the endpoint region one can extract the shape 
function and with it some fundamental matrix elements of higher-dimension 
operators in QCD. 

Starting from a resummation of the theoretical lepton spectrum, we pre- 
sent in Sec. ^ a heuristic argument that leads to the notion of a function 
S{y), which describes the fall-off of the spectrum in the endpoint region. In 
Sec. 0, we introduce this shape function, discuss its properties, and derive 
relations for the first two moments of S{y). Sec. ^ is devoted to a formal 
definition of the shape function to all orders in 1/mb. The leading contri- 
butions to the moments are related to forward scattering matrix elements of 
local, higher- dimension operators in HQET. For the purpose of illustration, 
a simple model calculation of the shape function is presented in Sec. ^. In 
Sec. H, we summarize our results, indicate possible further applications and 
improvements of the method, and give some conclusions. 

2 Resummation of the singular terms 

To motivate the concept of a shape function, let us try to resum the singular 
contributions in (|I|) into a corrected parton model decay distribution. Ob- 
viously, the term proportional to 5(1 — y) can be absorbed by a shift of the 
argument of the step function in the leading-order term. More interesting is 
the contribution proportional to 6'{1 — y). It arises at second order in the 
expansion of the step function. However, since there is no 5-function contri- 
bution of first order in 1/mf,, it follows that one needs more than one step 
function, resulting in a dispersion of the spectrum. In fact, to order I /ml, 
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we can rewrite the theoretical spectrum in the following way: 



N 



^ yF(y)l5:e(i-y + e,), (6) 

i=l 



2Tb dy ' N 



where 



r 1 \^ + 33A2 2 1 2 -^1 /^7^ 



h 



From the second relation it follows that the displacements Si are of order 
l/mfe. Thus, the first relation corresponds to a nontrivial cancellation. Note 
that ay can be identified with the characteristic width of the endpoint region, 
i.e., the region over which the spectrum is dominated by nonperturbative 
effects. This is already a remarkable conclusion: The coefficient of the most 
singular term in (|1|), which had no effect in the approaches of Refs. ^j, 
determines the size of the endpoint region. 

At this point, it is worthwhile to obtain some estimates of the nonper- 
turbative corrections. From the observed value of the mass splitting between 
B and B* mesons one obtains A2 — 0.12 GeV^. The parameter Ai, on the 
other hand, is not directly related to an observable. Recently, we have shown 
that the field-theory analog of the virial theorem relates the kinetic energy 
of a heavy quark inside a hadron (and thus Ai) to a matrix element of the 
gluon field strength tensor [|l6l. This theorem makes explicit an intrinsic 



"smallness" of Ai, which was not taken into account in existing QCD sum 
rule calculations of this parameter [jl^, [18], . As a consequence, we expect 
that (— Ai) is considerably smaller than predicted in these analyses. Here we 
shall use the range — Ai = 0.05 — 0.30 GeV^. According to its definition, Ai 
is negative, so that the width ay in (|^) is well defined. Using these numbers, 
as well as = 4.8 GeV, we estimate 5y ~ —0.03 and ay ^ 0.03 — 0.07. 
We can multiply these quantities by mb/2 to obtain the corresponding shift 
and spread in the lepton energy spectrum. They are 5E ~ —65 MeV and 
aE = (— Ai/12)^/^ ~ 65 — 160 MeV. The value of aE can be compared to 
the width of the gap between the parton model endpoint of the spectrum 
and the physical endpoint, which, if we neglect the pion mass, is located 
at EP^^'^ = mB/2. This width is AE ~ (mg - mfe)/2 ~ 240 MeV. These 
numbers seem quite reasonable. In fact, assuming that the distribution of 
the displacements Ei around y = 1 is approximately symmetric, we have to 
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require that < AE, which is equivalent to — Ai < 3{mB — mi,)^. For 
reasonable values of Ai, this bound is always satisfied. 

In Fig. |I|, we show the resummed lepton energy spectrum (^) for Ai = 
—0.2 GeV^, A2 = 0.12 GeV^, = 10, and a particular set of Ei satisfying 
the constraints in (^. For this choice of parameters, the dispersion in the 
spectrum is such that the endpoint falls close to the physical endpoint at 
y ~ 1.1. Our reinterpretation of the QCD-predicted lepton spectrum has led 
to a reasonable shape which, however, is quite arbitrary. In fact, increasing 
A^, we can generate any decreasing function that satisfies the constraints in 
(|^. In the following section, we will introduce a shape function S{y) instead 
of the sum over step functions. The constraints will then turn into predictions 
for the first two moments of this function. 



3 The shape function 

We proceed by replacing the sum over step functions in by a continu- 
ous function which we furthermore decompose as -d^y) = 9(1 — y) + 
S{y) F{l)/F{y). The form of the second term is chosen for later convenience. 
We shall refer to S{y) as the shape function. The support of this function 
is restricted to a small interval of width 2 A around y = 1, where A of or- 
der l/mb.0 Some basic properties of S{y) can be derived from the physical 
requirements that the differential decay rate by positive, and that i){y) be a 
continuous function. We note that 



S{y) = if |y-l|>A, 

Siy) > if y>l, 

lim ^(1 + e) - 5(1- e) = 1, 

limS'(l + e) - 5'(l-e) =0. (8) 

Moreover, we expect that 5"(y) < if y 7^ 1, but we shall not impose this as 
a condition on S{y). 

As emphasized in Refs. ^, because of the singular form of the operator 
product expansion one has to integrate the theoretical lepton spectrum with 

^More precisely, we should require that S{y) be exponentially small outside an interval 
of width 2 A. 
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a smooth function before it can be compared to data. On the set of smooth 
functions (i.e., functions of y which are slowly varying over scales of order 
1/mf,), a rapidly varying function such as Siy), which vanishes outside a 
small interval around y = 1 , obeys a singular expansion of the formQ 

^(y) = f: ^^^"Hi-y), (9) 

where the moments Mn are defined as 

oo A 

M„ = jdy (y - 1)" S{y) = Jdy{y- 1)" S{y) . (10) 

-A 

To see that (^ is correct, assume that any reasonable test function f{y) can 
be Taylor-expanded around y = I, with the result that 



/dy/(y)5(y) = E -f /^"Hl)- (11) 



In terms of the shape function, the theoretical lepton spectrum takes the 
form 

1 dr 



2Th dy ^ 



F{y)&{y) + F{l)S{y) . (12) 



From a comparison with (|T]), we find that the first two moments of S{y) must 
satisfy^ 



Notice that M„ = 0{l/'rn^'^^) by dimensional analysis. Hence, the QCD 
prediction that Mq = 0{l/m^) corresponds to a nontrivial cancellation: The 

■^This procedure is familiar from the multipole expansion of a localized distribution of 
charges in electrodynamics. 

^Our judicious choice of the prefactor y in front of S{y) has eliminated the contribution 
of Ai to Mq. 
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area under the shape function (almost) vanishes. The first moment, which is 
related to the characteristic width of the endpoint region, is of the expected 
order of magnitude. 

The concept of a shape function exploits to full extent the physical in- 
formation contained in the coefficients of the singular terms in the QCD- 
predicted lepton spectrum. We find that over a region of width 2 A oc l/m^, 
the spectrum is of genuinely nonperturbative nature and described by a func- 
tion S{y), the moments of which can be addressed in QCD. When one goes 
to higher orders in the l/rrib expansion, one can address higher moments. In 
fact, the moments obey an expansion of the form 

Af„ = ;^ + ;^ + .... (14) 

where so far we know the coefficients = 0, 6o = — IIA2/2, and ai = — Ai/6. 
With the exception of the moment Mq, where the leading term vanishes, 
we may argue that it would be a good approximation to know the leading 
coefficient a„ for each moment. The corrections involving bn only change 
the moments by small amounts. On the other hand, knowledge of a new 
moment teaches us a new piece of essential information about the shape 
of the spectrum in the endpoint region. The higher moments give a small 
contribution to integrated quantities such as the total decay rate, simply 
because in (|l^) one integrates over a small region. Nevertheless, they can 
affect the shape of the endpoint region in a substantial way. What is relevant 
to the shape are the rescaled moments 



h 

Mn= I dE, {2Ee - SiEe) = m^^ M„ = a„ + — + ... , (15) 







which remain nonzero in the limit rrif, — > 00. As an illustration of the im- 
portance of higher moments, we show in Fig. ^ two shape functions which 
have the same first two moments Mq and Mi, but different third (and higher) 
moments. The total decay rate and the average lepton energy are the same 
in both cases (up to terms of order 1/m^), but obviously the behavior close 
to the endpoint is quite different. 
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4 Formal definition of the shape function 



The above discussion shows that for an understanding of the lepton spectrum 
in the endpoint region, it is insufficient to truncate the theoretical calculation 
at order Instead, what one needs to investigate to all orders in 1/mb 

are the coefficients a„ in (JT^) and (p!5[) . They arise from the most singular 
terms in the shape function. In this section, we give a formal definition 
of these terms to all orders in l/m^. This will provide us with a relation 
between the coefficients a„ and forward scattering matrix elements of local, 
higher-dimension operators in HQET. 

As mentioned in the introduction, the derivation of the lepton spectrum 
is based on the operator product expansion in connection with an expansion 
of hadronic matrix element in powers of l/rrib, as provided by HQET. This 
is explained in detail in Refs. |^, |, ||, [l^l- Using the same technology, we can 



derive a closed expression for the most singular terms of the shape function, 
where S{y) is defined as in (|T2|) as the sum of all terms in the theoretical 



spectrum that become singular in the limit y ^ 1. We obtain the formal 
result 



s{y) = 



2 

1 — (v — p) -iD 

rrib 



6(1 — ?/))+ less singular terms, (16) 



which is valid to all orders in the l/rrib expansion. Here 

P = Pe/iTT'b denotes the rescaled lepton momentum, and we define the 
expectation value of an operator O as 

_ {Biv)\KOK\Biv)) 
^ ' {B{v)\KK\B{v)) ■ ^ > 

Here, hy is the velocity-dependent heavy quark field in HQET 0, and the 
states are the eigenstates of the corresponding effective Lagrangian. Details 
of the derivation of (0), as well as the extension to the case oi B ^ X^i-i' 
decays, will be given elsewhere [ pO|] . 

Expanding our result in powers of l/rrib, we obtain 

OO 1 / 2 \ " 

+ less singular terms. (18) 
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The forward scattering matrix elements between B mesons, or between any 
other hadronic states that are unpolarized, can be parameterized in the form 



( iD^^ . . . iD^") = AnV^^ ... w''" + terms with 51^'^^ (19) 

Since {v — pY = {1 — y) vanishes at the endpoint, only the coefficients An 
contribute to the most singular terms in S{y). For the same reason, we can 
replace factors 2v ■ {v — p) = {2 — y), which arise upon contraction of the 
indices in (|18D, by 1. This leads to the following expression for the shape 
function: 

1 A 

^iy) = — r — 5^""^^(1 - y) + less singular terms. (20) 



From a comparison with and (H), we obtain for the moments of S{y): 

(n + 1) m^^ n + 1 

The ffist three coefficients An are given by Aq = 1, Ai = 0, and A2 = — Ai/3. 



5 A model calculation 

It is instructive to consider a simple model for the shape function S{y). 
For this purpose, we evaluate the expectation value in (|16D following the 
phenomenological approach of Altarelli et al. (ACM) [Q. We emphasize, 
however, that this is mainly meant as an illustrative example rather than a 
prediction of the physical shape function. In fact, we will see very clearly the 
limitations and shortcomings of the ACM approach. 

In the ACM model, one assumes the validity of the parton model and 
incorporates bound state effect by assigning a momentum distribution 0(|pf,|) 
to the heavy quark inside the B meson at rest. It is then appropriate to 
replace the covariant derivative in (|1^) by the spatial components of the 
heavy quark momentum pb. The gluon field in the covariant derivative is 
neglected. Accordingly, in the rest frame of the B meson, one makes the 
replacement 

A(,_^).,zp^^ = ^eos^, (22) 

nib TTT'b 
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where is the angle between the lepton and the heavy quark momentum. 
Since we are interested in the behavior in the endpoint region, we can set 
y = 1. The matrix element in (|I^) is now replaced by an integral over the 
momentum distribution of the heavy quark: 



SACMiy) = /d|p,||p,p0(|p,|) |^^{e[i-i/ + ^cos^] -e(i-y)}. 

-1 

(23) 

It is straightforward to calculate the moments of this model shape function, 
and from (|2l|) the corresponding predictions for the hadronic matrix elements 
An. We find that A^^m _ 



oo 



A^^"" = { iPbi'i = J m\ bir^-^^) <pm) . (24) 



In the ACM model, one assumes a gaussian distribution. 



0(b1,|) = ^i^exp(-^^), (25) 



where pp is the Fermi momentum. This leads to the shape function 
SAcmiy) 



l-eil-y)]^(^\y-l\). (26) 



2 



^{x) = -y^ Jdte *^ denotes the complementary error function. For the coef- 
ficients A2n, one obtains 

<^" = ^^^^pr- (27) 

Comparing this to the general relation A2 = — Ai/3, we derive 

- A^c^ = IpI^ 0.08 GeV^, (28) 

where we have used — 230 MeV as obtained from the most recent fit of 
the ACM model to experimental data The model thus predicts a rather 
small value of — Ai. 
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In Fig. ^, we show the shape function of the ACM model for three different 
values of the Fermi momentum. We stress again that this simple model cal- 
culation is presented for pedagogical purposes only. In particular, note that 
in the ACM model the even moments of the shape function (corresponding 
to odd coefficients A2n+i) vanish by rotational invariance. This is a conse- 
quence of the fact that one replaces the operator of the covariant derivative 
by a c-number momentum vector. Whereas in QCD the commutator of two 
covariant derivatives gives the gluon field strength tensor, in the ACM model 
this commutator vanishes. However, exactly those terms involving the gluon 
field are responsible for an asymmetry in the shape function around y = 1. 
To see this, consider the matrix element of three covariant derivatives. Using 
the equations of motion of HQET, it is easy to show that 

{iD^'iD'^iD'^) = Asiv'^v'^ - g'"^)v\ (29) 

and taking the antisymmetric combination in /i and u, we find that is 
related to a matrix element involving the gluon field strength tensor: 

( ig, iD"" ) = A3 {v^'g'"' - v^'g^"') . (30) 

In QCD, there is no reason why such a matrix element should vanish. Hence, 
we expect an asymmetry of the physical shape function. The above example 
is instructive since it shows that a measurement of the moments of the shape 
function can provide quite fundamental information about the dynamical 
properties of the theory of strong interactions. 

6 Summary and conclusions 

We have presented a QCD-based approach to the inclusive lepton energy 
spectrum m B X„ I v decays. We have introduced the concept of a shape 
function, which is a genuinely nonperturbative object that describes the rapid 
fall-off of the spectrum in the endpoint region. The moments of this function 
obey a very simple relation to forward scattering matrix elements of local, 
higher- dimension operators. QCD predicts that the leading contribution to 
the first moment vanishes. The second moment, which is a measure of the 
size of the endpoint region, is proportional to the expectation value of the 
kinetic energy of the heavy quark inside the hadron. 
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Our approach goes beyond previous work on the lepton spectrum ||, |], 
which was apphcable only for lepton energies not too close to the endpoint. 
It aims at a systematic use of the rich source of information contained in the 
endpoint region. Whereas the main part of the spectrum is determined by 
kinematics and only receives small nonperturbative corrections, the behavior 
close to the endpoint is characterized by an infinite set of hadronic matrix 
elements. It is worth noting that, although there are no first-order power 
corrections to the spectrum and decay rate at small lepton energies, there 
exists a small region where the true spectrum is very different from the 
theoretical prediction. 

There are obvious improvements of the analysis presented here. Before 
confronting our results with data, it is necessary to include radiative correc- 
tions. For the case of -B — C v decays, they are known to affect the parton 
model spectrum in a significant way W\ Such corrections will affect the 



form of the shape function, too. We expect small perturbative corrections 
of order ^^(mf,) to all moments of the shape function. Moreover, radiative 
corrections will wash out the step in the shape function at y = 1, resulting 
in a rapidly varying, but not discontinuous, behavior. Another important 
generalization of our approach is that to the case oiB —^X^iv decays. The 
nonvanishing mass of the charm quark will lead to technical complications, 
but conceptually there is no problem in defining a shape function S'(?/, p) for 
any value of p = m^/m^. The (appropriately defined) moments of this gener- 
alized shape function are still related to the same hadronic matrix elements 
as in the case of a massless final state quark. Work on all these issues is in 
progress and will be reported elsewhere ||2^ . 

We believe that our approach will eventually lead to a better understand- 
ing of the nonperturbative aspects of inclusive decay spectra, with implica- 
tions for the measurement of some of the fundamental parameters of the 
standard model, such as the heavy quark masses and the elements Vch and 
Vuh of the Kobayashi-Maskawa matrix. Currently, the most promising appli- 
cations of the method seem to be the following: 

In B decays into charmless final states, an understanding of the end- 
point region is crucial for a reliable determination of Vuh- The approach of 
Refs. 0, H, 1^ cannot be used for this purpose, since current measurements 
are limited to a small energy range 2.3 < Ef < 2.6 GeV (corresponding to 
0.96 < y < 1.08 for = 4.8 GeV), which is too close to the endpoint. 
What is needed is some insight into the nonperturbative effects relevant to 



13 



the shape of the spectrum in the endpoint region. An expansion in powers 
of 1/mfc is not suitable for such a situation. The relevant physics is encoded 
in the moments of the shape function, which are related to forward scat- 
tering matrix elements of local, higher-dimension operators. Such matrix 
elements can be addressed using nonperturbative techniques such as lattice 
gauge theory or QCD sum rules. It may even be possible in these approaches 
to attempt a direct calculation of the shape function from its definition in 

For B Xciv transitions, the situation is very different. Already, there 
exist very accurate measurements of the lepton spectrum in this case. For a 
determination of Vcb and of the quark masses mi, and mc, an understanding 
of the endpoint region is thus not a necessary requirement. However, these 
decays offer the exciting possibility to extract the shape function from the 
data, simply by subtracting the (corrected) parton model spectrum from the 
measured distribution. One could then compute the first few moments of 
the shape function and extract some of the coefficients An-, which encode 
fundamental dynamical properties of QCD. The fact that QCD predicts the 
size of the first moment Mq in (|T3|) provides an important constraint, which 
can help to obtain a very precise determination of the 6-quark mass. 

We conclude with a speculation about yet another exciting possibility, 
namely to combine the analyses oi B — >■ XcCi' and B — >■ X^tV decays. 
In fact, one can derive a simple relation between the shape functions for 



these two processes [Q. One can imagine to measure the shape function 
in B decays into charmed particles, and then predict the shape function for 
charmless transitions. This avenue may well be the most promising one with 
respect to a precise extraction of Vuh- 
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Figure Captions 



Figure 1: An example of a resummed lepton energy spectrum according to 
(^. On the vertical axis, we show dT/dy in units of 2rf,. 



Figure 2: Two shape functions with identical moments Mq and Mi (a), and 
the corresponding lepton spectra (b). 



Figure 3: The shape function of the ACM model (a) and the corresponding 
lepton spectrum (b) for — Ai = 0.05 GeV^ (dashed), 0.1 GeV^ (solid), and 
0.2 GeV^ (dotted). The corresponding values of the Fermi momentum are 
Pp ~ 180 MeV, 260 MeV, and 365 MeV, respectively. In (b), the parton 
model spectrum is shown as a grey line. 
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